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Abstract — The chemical compositions of magmas can be derived from the compositions of their 
equilibrium minerals through mineral /magma partition coefficients. This method cannot be applied 
safely to basaltic rocks, either solidified lavas or cumulates, which have chemically equilibrated or 
partially equilibrated at subsolidus temperatures, i.e., in the absence of magma. Applying mineral/ 
melt partition coefficients to mineral compositions from such rocks will typically yield “magma com- 
positions” that are strongly fractionated and unreasonably enriched in incompatible elements (e.g., 
REEs). In the absence of magma, incompatible elements must go somewhere; they are forced into 
minerals (e.g., pyroxenes, plagioclase ) at abundance levels far beyond those established during normal 
mineral /magma equilibria. Further, using mineral /magma partition coefficients with such rocks may 
suggest that different minerals equilibrated with different magmas, and the fractionation sequence of 
those melts (i.e., enrichment in incompatible elements ) may not be consistent with independent con- 
straints on the order of crystallization. Subsolidus equilibration is a reasonable cause for incompatible- 
element-enriched minerals in some eucrites, diogenites, and martian meteorites and offers a simple 
alternative to petrogenetic schemes involving highly fractionated magmas or magma infiltration meta- 
somatism. 


1. INTRODUCTION 

The compositions of basaltic magmas are among the best 
available constraints on the compositions of their parent plan- 
ets and on petrogenetic processes within them (Drake, 1980). 
For bodies other than the Earth and Moon, basaltic meteorites 
are the only available basalt samples; however, meteorites are 
rarely large, commonly do not have bulk compositions rep- 
resentative of magma compositions, and commonly are brec- 
ciated into nonrepresentative fragments. Similarly, cumulate 
igneous rocks from all planets do not retain the compositions 
of their parent magmas. So there has been a strong impetus 
to derive the compositions of basaltic magmas from the com- 
positions of mineral grains that crystallized from them. In the 
main, these efforts have been successful because of a strong 
database of experimental studies on the distributions of ele- 
ments between magmas and minerals (e.g., Jones, 1995), and 
the development of extraordinarily sensitive analytical meth- 
ods, e.g., secondary ion mass spectrometry (SIMS) and syn- 
chrotron X-ray fluorescence (SXRF). 

A serious complication between a mineral analysis and the 
inference of its parental magma’s composition is postmag- 
matic (subsolidus or metamorphic ) redistribution of elements 
(Barnes, 1986a). A failure to recognize this process can lead 
to wildly erroneous estimates of trace element abundances in 
parent magmas, and thereby to completely unrealistic scenar- 
ios of magmatic and planetary evolution. The importance of 
subsolidus equilibration was first recognized by Consolmagno 
and Drake ( 1977), . . re-equilibration of the REE among 

phases means that attempts to calculate REE abundances in 
equilibrium with a single . . . phase will lead to erroneous 
conclusions . . .”. Subsolidus equilibration has been recog- 
nized in many relevant systems (e.g., Kennedy et al., 1992; 
Mittlefehldt et al., 1992; Fowler et al., 1994), but the char- 
acter of the erroneous conclusions has not been explored. 


In this work, I consider the effects of subsolidus equilibra- 
tion on magma compositions inferred from mineral analyses, 
emphasizing rare earth element (REEs) distributions in sim- 
ple basaltic systems. The derived principles are, however, rel- 
evant to all elements in all igneous rocks. First, I will explore 
the consequences of subsolidus equilibration on calculating 
magma compositions from a solidified basalt magma, then 
consider a basaltic cumulate, and finally reinterpret a few da- 
tasets from the literature. 

2. THEORY 

The partitioning of a trace element Z between two phases, 
1 and 2, is commonly described by the equilibrium distribu- 
tion coefficient D, the ratio of the mass concentrations of the 
element in the two phases at chemical equilibrium, 

Dz 2 = CyJCy. (1) 

Nomenclature here follows Beattie et al. ( 1993). If the par- 
tition coefficient and the element concentration in one phase 
are known, then the element concentration in the other phase 
can be calculated. For instance, if 1 were pigeonite pyroxene 
and 2 were basalt magma, knowledge of /)P'£ e,,n,tt ' hasd,, would 
allow one to calculate the concentration of Z in the basalt from 
the concentration in its equilibrium pigeonite. More generally, 
consider a system of unit mass composed of phases, 1,2, 
n , in mass proportions M {y M 2 , • * * A/ n . For a trace element 
Z with total system abundance ( = concentration) C/'\ Eqn. 
1 and mass balance on Z yield: 

n 

C?" = C/ - (A/ 1 + X M* 'Dz'). (2) 
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If C/ ot and all D and M values are known, C' A can be calculated 
directly, and the concentrations of Z in other phases can be 
calculated via Eqn. 1 . 

These equations are most commonly applied to igneous 
systems; if one phase is a magma, Eqn. 2 describes the dis- 
tribution of elements in batch melting or equilibrium crystal- 
lization. For example, consider the distribution of La, a light 
(highly incompatible ) rare earth element ( LREE) in a eucrite 
basalt containing C™ = 10 X CL If the basalt contains 50% 
magma, 25% pigeonite, and 25% plagioclase at chemical 
equilibrium, Eqn. 2 implies that C™ gma = 19.6 X Cl, 
C[”f com,c = 0 039 x c , and = 0 82 x CI using the 

partition coefficients of Table 1 . Magma is the main repository 
for highly incompatible elements like La. 

2.1. Subsolidus Equilibration 

In a basaltic rock that has chemically equilibrated at sub- 
solidus temperatures (i.e., in the absence of magma), all el- 
ements must be sited in solid mineral phases. Elements that 
are not abundant enough to form separate minor phases (e.g., 
the REEs) will become distributed among the solid phases. 
The REEs are concentrated in phosphates relative to basalt 
magma (Table 1), and similarly relative to solid pyroxene 
and plagioclase (Tables 2, 3). So the REE contents of indi- 
vidual phases in an equilibrated basalt depends on the bulk 
REE content and the mineral proportions. Other incompatible 
elements, like Zr, may form a separate minor mineral, like 
zircon or baddeleyite. In this case, abundances of Zr in other 
minerals may reflect saturation in the minor mineral and be 
independent of the bulk Zr content (i.e., the original magma 
composition ) . 

To calculate the distribution of elements in a basaltic rock 
that has equilibrated subsolidus, one may apply Eqn. 2, which 
was constructed without prescience of the phases involved. 


Table 1 . REE Mineral/Basalt Partition Coefficients 
(D), see Equation 1. 



pig/bas 

pig/bas 

aug/bas 

pl/bas 

wh/bas 

La 

0.002 

0.001 

0.03 

0.042 

25 

Ce 

0.004 

0.004 

0.04 

0.03 

25 

Pr 

0.012 

- 

— 

0.025 

25 

Nd 

0.019 

0.010 

0.10 

0.02 

25 

Sm 

0.031 

0.042 

0.17 

0.0175 

25 

Eu 

0.011 

— 

0.11 

1,2 

23 

Gd 

0.043 

— 

0.2 

0.016 

25 

Tb 

0.055 

- 

0.22 

0.0095 

25 

Dy 

0.07 

0.070 

0.24 

0.0085 

22 

Ho 

0.085 

- 

025 

0.008 

19 

Er 

0.1 

0.104 

0.27 

0.0075 

15 

Tm 

0.115 

- 

0.28 

0.007 

11 

Yb 

0.13 

0.129 

0.29 

0.0065 

8 

Lu 

0.13 

— 

0.30 

0.0068 

8 

Ref. 

1 

2 

1 

3 

4 


Basalt, bas; pigeonite, pig; plagioclase; augite, aug; 
whitlockite, wh. 

1 McKay et al. (1986a). 

2 Pun and Papike (1994). 

3 From and extrapolated from Phinney and Morrison (1990). 

4 From and extrapolated from McKay et al. (1986b, 1987) 
and McKay (priv. comm.) 


Table 2: REE Mineral/Mineral Partition 
Coefficients (D) for Magmatic Temperatures. 



pi/pig 

pl/aug 

pig/aug 

wh/aug 

La 

21 

1.4 

0.067 

830 

Ce 

7.5 

0.75 

0.10 

625 

Pr 

2.1 

— 

— 

— 

Nd 

1.0 

0.20 

0.19 

250 

Sm 

0.56 

0.10 

0.18 

150 

Eu 

110 

11 

0.10 

210 

Gd 

0.37 

0.08 

0.22 

125 

Tb 

0.17 

0.040 

0.25 

113 

Dy 

0.12 

0.035 

0.29 

92 

Ho 

0.09 

0.032 

0.34 

76 

Er 

0.075 

0.028 

0.37 

56 

Tm 

0.06 

0.025 

0.41 

40 

Yb 

0.05 

0.022 

0.45 

28 

Lu 

0.05 

0.022 

0.43 

27 


Calculated as quotients of values in Table 
1. Mineral abbreviations as in Table 1. 


To use Eqn. 2 for rocks equilibrated in the subsolidus, one 
needs mineral /mineral partition coefficients ( Tables 2, 3 ) . For 
magmatic and near-solidus temperatures, mineral /mineral 
partition coefficients can be estimated from mineral /magma 
coefficients. For minerals 1 and 2: 

rj I /magma 

Z 2 /magma * V ^ / 

(see Eqn. 1 ), provided the mineral /magma coefficients have 
been determined for identical or similar chemical and physical 
conditions and are applied to systems with similar mineral 
compositions. Partition coefficients for lower temperatures 
can be calculated from analyses of coexisting phases in met- 
amorphic rocks. Spot analyses (e.g., by focused beam instru- 
ments) are preferable, as whole mineral grains commonly 
contain inclusions of trace-element -rich materials. Electron 
microprobe analyses are abundant, but access only minor and 
the most abundant trace elements. Spot analyses by secondary 
ion mass spectrometry (SIMS) and synchrotron X-ray fluo- 
rescence (SXRF) are relatively rare at this time, but can in- 
clude the REEs. Analyses of mineral separates by neutron 
activation may include a wide suite of trace elements, but may 
be compromised by inclusions of other minerals (e.g., Bacon, 
1989). 

Mineral /mineral partition coefficients for magmatic con- 
ditions cannot be safely applied to subsolidus equilibria, be- 
cause of differences in temperature and because of phase 
transformations. The dependencies of mineral /mineral Ds on 
temperature have been explored in only a few systems, but D 
values tend to become more extreme with decreasing tem- 
perature (Table 3; Kurat et al., 1980). For instance, 

Df ugUe/orthopyroxc n e incfeases from ^45 at _ 115( )° C tO ~200 

at '-~1050°C in terrestrial mantle xenoliths, based on mineral 
separates ( Stosch, 1982). 

D values also vary with mineral composition (e.g., Coulson 
et al., 1988, 1989; Gallahan and Nielsen, 1992) and this vari- 
ation undoubtedly extends to subsolidus equilibration. The 
effect of mineral composition on mineral /mineral Ds is most 
significant for pyroxenes, particularly the effect of their cal- 
cium content. Because Ca “opens” the pyroxene structure to 
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permit the larger trace cations to enter (e.g., McKay et al., 
1986a; Shearer et al., 1989), more calcic pyroxenes will tend 
to contain more of the REEs, especially LREEs (and Al, Ti, 
etc.) . Thus, abundances of larger trace cations in pyroxenes 
will tend to be ordered 

(Shearer et al., 1989), implying that DtT c,p '^ mixc 
< OtT' L/orlhopyroxcnc ■ This tendency is not perfect; D ^ omlc/dU £' lc 
is ~0.4 under magmatic conditions (Table 2), and so is 
^(mhc.pymxcne/augiie j n man y metamorphically equilibrated eu- 
crites and lunar samples (Table 3; Floss, 1991; Pun and 
Papike, 1995). Other compositional variables (e.g., Fe/Mg, 
Al and Ti content of pyroxene) may also affect D values 
(Coulson et al., 1988, 1989; Gallahan and Nielsen, 1992). 

2.LL Example: Near-solidus equilibration of basalt 

Consider again a basalt rock (e.g., a eucrite) with CS, 
= 10 X Cl (Fig. 1 ). If it equilibrated at a subsolidus but near 
solidus temperature, it might consist of 49.5% plagioclase, 
49.5% pigeonite, 0.1% whitlockite (a Ca phosphate), and 
0.9% minor phases with negligible REE abundances (e.g., 
silica, Fe metal, troilite); the exact mineral proportions are 
not critical. Figure 1 shows the REE contents of the basalt 
minerals, as calculated using the partition coefficients of Ta- 
ble 1 . The REE abundances in these pigeonite and plagioclase 
samples are much higher than those in equilibrium with basalt 
magma. In the example above of a basalt (CJkh — 10 X Cl) 
with 50% magma, one calculates C£f comlc = 0.04 X Cl and 
^ piagiociasc = Q g x CI; with no magma present, C[“ = 0.4 

X Cl and CJ’ l a agK,clasc = 9 X Cl. This tenfold enrichment in 
mineral REE abundances reflects the absence of magma. 
When magma is present, the magma is a major carrier of the 
REEs (and other incompatible elements); when magma is 
absent, the REEs must be accommodated by solid phases. 

Calculation of the hypothetical ‘ ‘parent basalt’ ’ by mineral/ 
melt partition coefficients ( Eqn. 1 ) is unjustified here, because 
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Fig. 1. Calculated REE abundances in minerals of a basalt equil- 
ibrated in the absence of magma just below its solidus temperature 
(D values of Table 2), and the hypothetical “parent magma” inferred 
by Dffi iaWb “ h from those mineral compositions (O values of Table 
1). See text for details of model. Whole-rock (heavy line) and equi- 
librium mineral compositions (Pig., pigeonite; Plag., plagioclase; 
Whit; whitlockite). 


the calculation assumes that the minerals last equilibrated in 
the presence of a basaltic magma. Following this false as- 
sumption (using D values of Table 1), one would infer a 
“parent basalt” with REE abundances as shown in Fig. 1. 
Because these same D values were used to calculate the REE 
abundances in the minerals and the hypothetical magmas, all 
minerals in this rock appear to have equilibrated with the same 
“parent basalt,” 

y'-'Tot 

C? ™""=- . (4) 

£ M, D 


Table 3. Mineral/mineral Partition Coefficients (D) for Subsolidus Temperatures. 



i pl/opx ! 

pl/aug 

opx/aug | 

wh/opx 

Sample 

76535 

62255 

62255 

62255 

60055 

Moama 

UsMo 

Roda 

La 

150 

7.5 

0.52 

0.07 

- 

0.010 

0.005 

500.000-1,500,000 

Ce 

30 

3.6 

0.23 

0.06 

0.006 

0.008 

0.005 

150,000-1,000,000 

Pr 

12 

17 

0.11 

0.07 

- 

- 

- 

- 

Nd 

5.3 

0.85 

0.06 

0.07 

0.011 

0.020 

0.0075 

-10,000 

Sm 

1.4 

0.31 

0.02 

0.08 

0.007 

0.042 

0.016 

-10,000 

Eu 

- 

- 

11* 

- 

- 

0.0006 

0.022f 

-85,000* 

Gd 

0.46 

0 13 

0.013 

0.11 

0.04 

- 

0.027 

- 

Tb 

0.20 

0.06 

0.009 

0.16 

0.07 

0.138 

0.036 

- 

Dy 

0.1 

0.03 

0.006 

0.19 

0.12 

- 

0.048 

- 

Ho 

0.07 

0.03 

0.008 

0.27 

0.14 

- 

0.059 

- 

Er 

0.04 

0.009 

0.003 

0.32 

0.23 

0.28 

- 

- 

Tm 

0.03 

- 

- 

0.38 

0.32 

- 

0.12 

- 

Yb 

0.04 

0.008 

0.004 

0.53 

0.43 

0.387 

0.14 

- 

Lu 

_ 

- 

- 

0.71 

0.46 

- 

0.19 

- 

Ref. 

1 » 

1 

1 

1 

2 

3 

4 


Mineral abbreviations as in Table 1, and: orthopyroxene, opx. 

1 Derived from data in Floss (1991). 

2 Pun and Papike (1995). 

3 Derived from data in Stosch (1982). 

4 Derived from data in Mittlefehldt (1994a) and Fowler et al. (1995). 

* Equilibration under reducing conditions, oxygen fugacity approximately at iron-wtistite. 
t Equilibration under oxidizing conditions, oxygen fugacity near fayalite-magnetite-quartz. 
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This hypothetical “parent basalt” (Fig. 1 ) has high overall 
abundances of the REE ( 130-200 X Cl), moderate enrich- 
ment in light REEs (La/Lu = 1.6 X Cl), and strong depletion 
in Eu, relative to Sm and Gd. 

Why is this hypothetical “parent basalt” (Fig. 1) so dif- 
ferent from the original simple basalt, with CJSU = 10 X Cl? 
The calculated Eu content of the “parent basalt” is similar to 
the original basalt, because Eu is partitioned into plagioclase 
in both subsolidus and super-solidus conditions. The remain- 
ing REEs, having no melt phase in which to reside, have been 
forced into solid phases; the rock contains so little whitlockite 
that the bulk of the REEs are in plagioclase and pigeonite. 
Thus, the plagioclase and pigeonite have been forced to ac- 
commodate far more of the REEs than they would have in the 
presence of magma. And a higher REE content in the plagio- 
clase and pigeonite require a higher REE content in the in- 
ferred hypothetical “parent basalt.” 

Thus, the geochemist or petrologist will be grossly misled 
by assuming or inferring that the compositions of minerals in 
this rock reflect equilibrium with magma. The inferred hy- 
pothetical “parent basalt” has — 20 times the REE content of 
the original basalt, is moderately fractionated while the orig- 
inal basalt is not, and has a strong Eu anomaly which the 
original basalt does not (Fig. 1 ). Certainly, the petrogenesis 
of this hypothetical “parent basalt” would be fundamentally 
different from the petrogenesis of the original simple basalt. 

2./. 2. Example : Deep-subsolidus equilibration of basalt 

If this basalt rock equilibrated at a lower temperature, sig- 
nificantly below the solidus, the results of applying mineral/ 
melt REE partitioning might be still more extreme, as sug- 
gested by a comparison of solidus-temperature vs. subsolidus 
D values (Tables 2,3). The effects of subsolidus equilibration 
on inferred magma compositions can be seen by again con- 
sidering a basalt with = 10 X Cl. If this basalt chemi- 
cally equilibrates at deep-subsolidus temperatures (e.g., 
700°C), it might consist of 49.5% plagioclase, 44.5% ortho- 
pyroxene, 5% augite, 0.1% whitlockite, and 0.9% other 
phases (again, the phase proportions are not critical). The 
augite and orthopyroxene derive from inversion of the origi- 
nal magmatic pigeonite. REE abundances in these minerals, 
shown in Fig. 2a, are calculated from equation 2 and the D RHK 
of Table 3. For consistency, I use the Ds for the lunar sample 
62255, except for which must be taken from 

Table 2 in the absence of an independent subsolidus estimate. 
Inverting those abundances to magma compositions (Z> values 
of Table 1 ), one finds that all the minerals would have formed 
from fractionated magmas much richer in REEs than the orig- 
inal bulk rock and magma (Fig. 2b). 

In addition, each mineral would appear to have equilibrated 
with a different magma composition (Fig. 2b) except for 
whitlockite and augite. Whitlockite and augite would appear 
to have come from the same magma because the calculation 
had to rely on f or magmatic temperatures. The 

relative REE enrichments in these hypothetical parent mag- 
mas derive entirely from the different temperature depend- 
encies of the D values, and have nothing to do with magmatic 
processes. 

For example, the order of increasing REE content of “par- 
ent magmas” could be wrongly construed as the order in 


which the minerals crystallized (in a fractionating system). 
Even though the plagioclase “parent magma” is poorer in 
REE than the pyroxene “parent magma” (Fig. 2b), plagio- 
clase and pyroxene co-crystal lize in the eucrite basalts on 
which this model was based. Looked at more quantitatively, 
the “parental magma” for plagioclase has roughly one-third 
the REE content of the “parent magma” for orthopyroxene. 
For closed system crystallization, this would imply crystalli- 
zation of approximately two-thirds of the system as plagio- 
clase before pyroxene began to crystallize. This inference 
must be wrong, as the system has only 49.5% plagioclase. 

2.1.3. Example: Subsolidus equilibration of a basaltic 
cumulate 

Cumulate igneous rocks do not ( ipso facto) retain the com- 
positions of their parent magmas, and so are natural targets 
for retrieving magma compositions by mineral /magma D val- 
ues. If the cumulate minerals have not chemically equilibrated 
at subsolidus temperatures, mineral /magma Ds are an effec- 
tive indispensable tool in recovering magma compositions. If 
the cumulate minerals have chemically equilibrated in the ab- 
sence of magma, application of mineral /magma D values will 
again suggest incorrect compositions for parental magmas. 

Consider as examples two nearly pure pigeonite orthocu- 
mulates. Let the first comprise 0.1% eucritic magma (REEs 
at 10 X Cl) and 99.9% pigeonite pyroxene that equilibrated 
with the magma; let the second comprise 1 % eucritic magma 
and 99% pigeonite. Let these cumulates crystallize and chem- 
ically equilibrate as closed systems at temperatures just be- 
neath the basalt solidus. After solidification, the first cumulate 
would contain —99.9495% pigeonite, —0.0495% plagioclase, 
—0.0001% whitlockite, and —0.0009% phases without sig- 
nificant REEs; the second cumulate would contain —99.495% 
pigeonite, —0.495% plagioclase, —0.001% whitlockite, and 
—0.009% phases without significant REEs. 

Figure 3 shows the REE contents of “magmas” for these 
two cases calculated from pigeonite compositions and 

D ^>ni,c/basal, Q f J J n bmh 

cases, 0.1% parent magma 
and 1 % parent magma, the HREE contents of the calculated 
“magmas” are within uncertainty of the true parent magma. 
But LREE abundances in the calculated “parent magmas” 
are significantly different. For the cumulate with 0.1% parent 
magma, the calculated “magma” appears to have La/Lu 
= 1.4 X Cl, which is only slightly above that of the true parent 
magma. Still, this elevated La/Lu would be sufficient to sup- 
port hypotheses of fractionated parent magma or of magmatic 
infiltration metasomatism. For the cumulate with 1% parent 
magma, the calculated “magma” appears to be highly en- 
riched in LREEs, with La/Lu = 4 X Cl. This strongly frac- 
tionated “magma” would imply a completely different pet- 
rogenesis than would the true eucritic parent magma. 

2.1.4. Conclusion 

Spectacular failure may attend the derivation of magma 
compositions from equilibrated rocks through application of 
mineral /melt partition coefficients. For basaltic systems, the 
derived magma compositions will be unnaturally enriched in 
incompatible elements and highly fractionated. Further, each 



REE /Cl -B REE /Cl 
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Fig. 2. Calculated REE abundances in minerals of a basalt equil 
ibrated at a temperature far below its solidus, and inferred hypothet- 
ical “parent magma” compositions from those minerals. See text for 
details of model, (a) REE abundances in a basalt and its constituent 
minerals (D values of Table 3, except for /?Khi- l,K ' kltc/aug,,c and 
£ > augiic/pige«>n,tc f rom T a t>] e 2). Whole rock (Bulk, heavy line) and equi- 
librium mineral compositions: Opx., orthopyroxene; Aug., augite; 
Plag., plagioclase; Whit, whitlockite. (b) REE abundances in hypo- 
thetical “parent magmas” inferred from mineral compositions of 
part a) and DS rill/hM " of Table 1. Whole rock (Bulk, heavy line) and 
hypothetical “magma” compositions in equilibrium with minerals: 
Opx. Magma, from orthopyroxene; Aug. + Whit. Magma, from au- 
gite and whitlockite; Plag. Magma, from plagioclase. 

mineral may yield a different magma composition, depending 
on the mineral species and the temperature of equilibration. 
The inference of different magma compositions during the 
crystallization of the rock may lead to spurious conclusions 
about order of crystallization and causes for changing magma 
composition. 

3. EXAMPLES FROM THE LITERATURE 

As shown above, a failure to recognize subsolidus or meta- 
morphic equilibration of trace elements can have dire geo- 
chemical consequences. A few cases from the literature will 
suffice as examples. 

3.1. Cumulate Eucrites 

Among the eucrite ( basaltic ) meteorites are a few cumulate 
lithologies, with relatively low REE abundances (—0.5-5 


X Cl), Eu excesses, and depletions of LREEs (e.g., Dodd, 
1981; Hamet et al., 1978). These characteristics suggest that 
the cumulate eucrites are plagioclase-pyroxene cumulates, but 
the REE characteristics of their parent magmas are in doubt 
(e.g., Consolmagno and Drake, 1977; Stolper, 1977; Hamet 
et al., 1978, Pun and Papike, 1995). 

Many investigators have assumed or inferred that mineral 
compositions in the cumulate eucrites reflect mineral /magma 
equilibria, and so have inferred or suggested that their parental 
magmas were highly fractionated (Ma et al., 1977; Stolper, 
1977; Ma and Schmitt, 1979; Pun and Papike, 1995; Hsu and 
Crozaz, 1995). However, the compositions of cumulate eu- 
crite minerals probably reflect mineral /mineral, not mineral/ 
magma, equilibria (Schnetzler and Philpotts, 1969; Consol- 
magno and Drake, 1977; Treiman and Drake, 1985), and re- 
cent SIMS data show that REEs in their pyroxenes were re- 
distributed at subsolidus temperatures (Pun and Papike, 
1995). 

For example, consider the Moore County cumulate eucrite. 
Figure 4 shows the REE contents of a Moore County plagio- 
clase separate, the estimated average pyroxene (Schnetzler 
and Philpotts, 1969; Pun and Papike, 1995), and of hypo- 
thetical magmas in equilibrium with both minerals. These hy- 
pothetical magmas are strongly enriched in REEs and espe- 
cially enriched in LREEs compared to normal eucrites (Fig. 
4), as shown by Ma et al. ( 1977), Stolper ( 1977), Ma and 
Schmitt ( 1979), and Pun and Papike ( 1995). However, the 
REE patterns of the minerals and inferred magmas in Fig. 4 
are qualitatively similar to those from equilibrated basaltic 
rock in Figs. 1 and 2. This suggests that the REE distribution 
between minerals reflects subsolidus equilibrium, as is con- 
sistent with the REE analyses and two-pyroxene thermometry 
of Pun and Papike ( 1 995 ) . The near-identity of the ‘ 'magma' ’ 
compositions inferred from plagioclase and pyroxene sug- 
gests that intergrain diffusion of REEs ceased at near-solidus 
temperatures, even though intragrain redistribution of REEs 



Fig. 3. Calculated REE abundances in two cumulates from a eu- 
critic basaltic parent magma (REE at 10 X Cl): 1 ) 99% pigeonite and 
1% parent magma; and 2) 99.9% pigeonite and 0.1% parent magma. 
For each case, figure shows calculated composition of bulk cumulate 
and composition of “parent magma” calculated from chemically 
equilibrated pigeonite (D Q f Table 2) in the cumulate and 
of Table 1 . 
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Ftc. 4. REE abundances in Moore County minerals and inferred 
“parent magmas” ( D values from Table 1). Plag., plagioclase, open 
circles (Schnetzler and Philpotts, 1969); Px, pyroxene* filled trian- 
gles* (weighted average from analyses of Pun and Papike, 1995). 


continued to much lower temperature (Pun and Papike, 
1995). 


3.2. Diogenites 

The diogenite meteorites are orthopyroxenites with minor 
chromite and olivine, and are commonly interpreted as igne- 
ous cumulate rocks associated in some way with the eucrite 
basalts (Dodd, 1981; Mittlefehldt, 1994a; but see Jurewicz et 
al., 1995). Almost all have nearly homogeneous mineral com- 
positions, and so are candidates for subsolidus chemical equil- 
ibration; Garland is the one exception, and Roda and Mane- 
gaon show considerable inter-grain chemical variations (Mit- 
tlefehldt, 1994a; Fowler et al., 1994, 1995). Within the 
candidate equilibrated diogenites, Fe/Mg is effectively con- 


stant, and there are limited ranges in Al, Cr, Ti, Zr, and REE 
contents (Mittlefehldt, 1994a; Fowler et al., 1994, 1995). 
Equally important, their orthopyroxenes have nearly constant 
Cr/Al abundance ratio (Fowler et al., 1994), which is not to 
be expected in a magmatic system crystallizing orthopyroxene 
± chromite (Barnes, 1986b; Coulson et al., 1988; Jurewicz et 
al., 1995). Thus, the homogeneity of Fe/Mg and apparent 
nonigneous element distributions in diogenites (except Gar- 
land) suggest that subsolidus chemical equilibria was important. 

3.2. L LREE-rich whitlockite 

In the Roda diogenite, Mittlefehldt ( 1994a) reported min- 
ute grains of REE-rich whitlockite with La = 65,000 X Cl 
( 1.8% La 2 0 1 ) and Sm = 6600 X Cl (0.1 1% Sm 2 O0 by elec- 
tron microprobe spot analysis (Fig. 5). If this whitlockite 
composition reflected equilibrium with a magma, it would 
have had La = 2600 X Cl and Sm = 260 X Cl (D values of 
Table 1 ) . A similar REE-enriched phosphate, merrillite, is pres- 
ent in the Mt. Padbury mesosiderite, although chemical equilib- 
rium among its phases is not complete (Kennedy et al., 1992). 

However, the REE-enrichment of the Roda whitlockite 
could be a result of subsolidus chemical equilibrium with or- 
thopyroxene. To test this interpretation, I will use 
D whi,UKki te / t mhopyroxene from subso i idus partition data, and com- 
pare those values to o^ kKkitcA,rthopymKene from Roda. 

There is considerable uncertainty in DREE OCkite/orthop yr oxcnc 
from subsolidus data, because there appear to be no indepen- 
dent measurements of this mineral /mineral distribution coef- 
ficient. DSf kl,c/urthopyroxene for subsolidus temperatures must 
be estimated as 

whulockite/tirthopyroxcnc _ ^jwhitUifkitc/augilCyi^jonhopyroxene/augite (5) 

using the range of D^ ,pyrmenc/aug,tc for subsolidus tempera- 
tures (Table 3) and the single value for DREE locki,c/auei,c for 



La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu 


Fig. 5. Whitlockite in the Roda diogenite. A comparison of D^h <x ' k,,c/or,hopyru<ent as measured in Roda and as predicted 
from the data of Tables 2 and 3 (see text). Dark-shaded field encloses permissible DKHK <xk,tc/onh<,pym, ' enc based on electron 
microprobe analysis of whitlockite in Roda and various SIMS and INAA analyses of Roda orthopyroxene (Mittlefehldt, 
1994a; Fowler et al., 1995). Light-shaded field encloses permissible f rom in Table 2 

and from Xable 3 
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magmatic temperatures (Table 2). The REEs are predicted to 
be very strongly partitioned into whitlockite; the range of D 
values is shown in Fig. 5. 

The natural values of o^ : f K ' k,,c/,m ^ ,pyroxerM: calculated from 
Roda must be limited to the light REEs (La, Ce, Nd, Sm, Eu ) 
because heavier REEs were not be detected by EMP (Mittle- 
fehldt, 1994a). Unfortunately, for Roda 

are uncertain because orthopyroxene adjacent to the whitlock- 
ite has not been analyzed, and because REE contents of Roda 
orthopyroxenes are quite variable (Fowler et al., 1995). Av- 
erage SIMS analyses for Ce shows a large range, 0.058 to 
0.34 X Cl (Fowler et al., 1995), and the other LREEs are 
also variable: La ranges from 0.047 X Cl (average SIMS spot 
analysis; Fowler et al., 1995) to 0.136 X Cl (INAA analysis; 
Mittlefehldt, 1994a). Sm only ranges from 0.58 X Cl (aver- 
age SIMS spot analysis; Fowler et al., 1995) to 0.64 X Cl 
(INAA analysis; Mittlefehldt, 1994a). The Mittlefehldt 
( 1994a) analyses, being on a chip sample, could have been 
compromised by inclusions of whitlockite. It is, of course, not 
known if the full range of orthopyroxene compositions have 
been sampled. 

As shown in Fig. 3, the predicted o^ ,,Kk,,c/ortbopyrt,xenc are 
in good agreement with the Roda for Sm, 

Nd, and Eu, but are an order of magnitude too low for La and 
Ce. Given the paucity of data and uncertainty about the REE 
content of orthopyroxene near the whitlockite, any number of 
hypotheses could explain the disagreement in Fig. 5 at La and 
Ce. Possibly, the whitlockite grew from a LREE-enriched 
“metasomatic” magma and never equilibrated with Roda py- 
roxenes (Wadhwa and Crozaz, 1994). Or, the orthopyroxene 
associated with the phosphate could be significantly richer in 
REEs than any Roda orthopyroxene analyzed so far. Or, the 
predicted d^***''«*w«'™* f or La and Ce could be too low, 
perhaps because the solidus-temperature £>RKH <Kk,,c/aueilc cannot 
be safely extrapolated to subsolidus temperatures. I favor the 
last explanation, as the first two seem inconsistent with agree- 
ment between observed and predicted o^ i ^ lc, ^ hopywxetK for 
Sm, Nd, and Eu. Thus, it seems reasonable to infer that REE 
abundances in Roda’s whitlockite represent chemical equilib- 
rium with adjacent orthopyroxene, and that Roda’s whitlock- 
ite and orthopyroxene can be used as a preliminary calibration 
for subsolidus (Table 3). 

3.2.2. REE-rich plagioclase 

Plagioclase is rare in diogenites, but a few large grains have 
been recovered from the Johnstown diogenite by Mittlefehldt 
( 1979, 1994a). These plagioclases are rich in REEs and par- 
ticularly the LREEs (e.g., La = Ce = 20 X Cl, La/Sm = 1.3 
X Cl). They were interpreted as having equilibrated with a 
magma that was highly enriched in the REEs, La = 270 X Cl, 
and strongly fractionated La/Sm = 3-7 x Cl (Mittlefehldt, 
1979, 1994a) . There is little other evidence for REE-enriched 
and fractionated magmas in diogenite petrogenesis, so Mittle- 
fehldt (1979, 1994a) suggested that the plagioclase grains 
might be exotic, and not related to the diogenite. 

The constraints of subsolidus equilibration are consistent 
with the plagioclase grains being exotic, but not necessarily 
with their origin in a LREE-enriched magma. An exotic origin 
is likely because the plagioclase composition cannot reason- 


ably result from equilibration with a diogenitic parent magma 
(Mittlefelhdt, 1979) nor from subsolidus equilibration within 
a diogenite. To show the latter point, consider a diogenite 
consisting of 98% orthopyroxene, 0.5% plagioclase, and 1 .5% 
other phases (chromite, troilite), with La = 0.33 X Cl and 
Sm = 0.40 X Cl (as does Johnstown; Mittlefehldt, 1994a). 
Using the plagioclase /orthopyroxene partition coefficients of 
Table 3, the equilibrium plagioclase would have La = 2.4- 
29 X Cl and La/Sm = 25-30 X Cl (ranges reflect use of 
coefficients from samples 76535 vs. 62255). The predicted 
La abundances are consistent with the Johnstown plagioclase, 
but the predicted La/Sm ratios are an order of magnitude too 
large. Impurities in the plagioclase sample will not help; py- 
roxene contains too little REEs, and whitlockite is likely to 
have La/Sm near 10 X Cl (see above; Mittlefehldt, 1994a). 

However, the Johnstown plagioclase need not have formed 
in a LREE-enriched and fractionated magma; it could reason- 
ably be a fragment from an equilibrated eucrite basalt. Eucritic 
fragments are present in other diogenites (Aioun el Atrouss, 
Garland, Pecklesheim: Delaney et al., 1983; Mittlefehldt, 
1994a), and are not unexpected, as the HED meteorite brec- 
cias include the full range of mixtures of diogenitic and eu- 
critic fragments (Delaney et al., 1983). The degree of frac- 
tionation of the Johnstown plagioclase (La/Sm = 1.3 X Cl) 
is similar to that of the calculated plagioclases in Figures 1 
and 2a (La/Sm = 2.5 and 2), but its absolute REE abun- 
dances are higher. The difference in absolute REE abundances 
could reflect a different temperature of equilibration, a parent 
rock more enriched in REEs than the model, or traces of phos- 
phate in the plagioclase. If the plagioclase sample of Fig. 2a 
contained —0.25% whitlockite, it would analyze at La = 20 
X Cl and La/Sm = 1.4 X Cl, nearly identical to that of the 
Johnstown plagioclase. 

3.3. Martian Meteorites 

Twelve of the basaltic achondrite meteorites are almost cer- 
tainly from Mars (McSween, 1994; Yanai, 1995; Mason, 
1995 ). These martian meteorites have become indispensable 
in the current understanding of Mars, especially in Mars’ ig- 
neous petrogenesis, and there is a great impetus to retrieve 
magma compositions from them (e.g., Stolper and McSween, 
1979; Longhi and Pan, 1989; Treiman, 1993). However, most 
of the martian meteorites are crystal cumulates, and most of 
these have experienced some degree of subsolidus chemical 
equilibration (e.g., McSween, 1985; Treiman et al., 1994; 
Mittlefehldt, 1994b). For these equilibrated and partially 
equilibrated meteorites, one might assume that minerals retain 
their magmatic compositions, and thereby possibly fall into 
error. 

3.3.1. ALH84001 

The ALH84001 meteorite is an orthopyroxenite, recently 
recognized by Mittlefehldt ( 1994b) as martian and not a di- 
ogenite. ALH84001 consists of approximately 97% orthopy- 
roxene, 2% chromite, 1% maskelynite, and 0.15% phosphate 
(Dreibus et al., 1994; Wadhwa and Crozaz, 1994). The min- 
erals in ALH84001 are essentially homogeneous in compo- 
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sition ( Mittlefehldt, 1994b; Papike et al., 1994), excepting 
those in the postigneous carbonate-rich alteration assemblage. 

Despite the chemical homogeneity of the minerals in 
ALH84001, the REE contents of parent magmas have been 
calculated from mineral /melt partition coefficients and SIMS 
analyses of minerals (Papike et al., 1994; Wadhwa and 
Crozaz; 1994). Orthopyroxene was inferred to have equili- 
brated with a magma that was slightly depleted in LREEs 
compared to HREEs (with HREEs ~10 x Cl, Papike et al., 
1994). Wadhwa and Crozaz ( 1994) further inferred that pla- 
gioclase and phosphate in ALH84001 equilibrated with mag- 
mas that were strongly enriched in REEs overall and espe- 
cially the LREEs, and postulated that ALH84001 was infil- 
trated during late crystallization by very highly fractionated 
magmas. 

Available evidence suggests, however, that REEs in 
ALH84001 minerals were affected by subsolidus equilibra- 
tion, and cannot be used safely to infer magma compositions 
or magmatic processes. As noted above, minerals in 
ALH84001 are effectively homogeneous. Further, it appears 
likely that REE distributions in the ALH84001 minerals re- 
flect subsolidus equilibria. From the limited data in Wadhwa 
and Crozaz ( 1994), one can calculate some REE distribution 
values for ALH84001 and compare them with the equilibrium 
values of Table 3. The data of Wadhwa and Crozaz ( 1994) 
imply that, for ALH84001, oS u,wkite/,,rthopyro,lcr * e « 10 4 5 and 
D whitiix.ki«e/onhop y roxtnc w 15; va i ues m comparable to 

(though slightly more extreme than) the equilibrium values 
of Tables 2 and 3 and Fig. 5. The Wadhwa and Crozaz ( 1994) 
analyses of maskelynite also imply £)^‘ (Klase/orth °pyro«nc w 75 
in the range of subsolidus equilibrium values of Table 3 (it is 
possible, however, that the analyzed maskelynite was contam- 
inated with shocked mesostasis material, which could be rich 
in REEs). Thus, one may credibly hypothesize that the ob- 
served REE distributions within ALH84001 reflect subsolidus 
equilibria, and have no significance for magma compositions 
or magmatic processes. 

4. AVOIDING PERIL: RECOGNIZING 
SUBSOLIDUS EQUILIBRATION 

As shown, applying mineral /magma partition coefficients 
to chemically equilibrated rocks can yield unrealistic magma 
compositions. So, how can one avoid this peril? How can one 
recognize rocks (or minerals) which might have experienced 
subsolidus equilibration? 

The most general warning sign of subsolidus chemical 
equilibration in an igneous rock is homogeneous mineral 
compositions, because extensive subsolidus equilibration will 
erase any chemical zoning imposed during igneous fraction- 
ation. However, homogeneous mineral compositions do not 
necessarily require subsolidus equilibration. Crystals may 
have grown in a magma chamber so large and well-stirred 
that their growth did not fractionate the magma significantly. 
Or, diffusion may be so rapid (e.g., Fe, Mg in olivine) or 
crystallization times so long (e.g., deep plutonic settings) that 
minerals may continuously equilibrate with a fractionating 
magma. In the latter case, however, one might expect subsol- 
idus cooling to be as protracted as the igneous cooling, and 
thus expect that subsolidus chemical effects might be signif- 
icant. 


Even if minerals are chemically zoned, they may still have 
been significantly affected by subsolidus chemical equilibra- 
tion. One might suspect subsolidus equilibration if the zoning 
is monotonic, and especially if the zoning is not consistent 
with igneous fractionations. Monotonic zoning itself (e.g., 
consistently increasing Fe/Mg from core to rim of a pyroxene) 
is not diagnostic, because it can develop during igneous frac- 
tionation as well as during metamorphic equilibration. How- 
ever, fine-scale oscillatory zoning does not develop during 
metamorphism, and can be taken as prima facie evidence that 
igneous compositions are retained (e.g., Treiman and Essene, 
1985). Second, any chemical zoning must be consistent with 
igneous fractionation patterns, i.e., mineral /magma distribu- 
tion coefficients. Consider, for instance, a basaltic magma 
crystallizing pyroxene and olivine. Chromium is compatible 
in pyroxene (D > 1 ) and Al is incompatible (D < 1 ) ; neither 
element enters olivine significantly. Thus, the Cr/Al ratio in 
the crystallizing pyroxene should decrease as pyroxene 
grows, as is seen in experimental studies of eucrite and di- 
ogenite genesis (e.g., Jurewicz et al., 1995 ) . In such a basaltic 
rock, pyroxenes with constant Cr/Al should be viewed with 
a suspicion of subsolidus equilibration. 

On the other hand, magma compositions inferred from min- 
eral/magma distribution coefficients can point to subsolidus 
chemical equilibration. Minerals affected by subsolidus 
chemical equilibration will tend to imply highly fractionated 
magma compositions, with high abundances of incompatible 
elements. Inference of such a fractionated evolved magma, in 
a setting where depleted or primitive magmas should prevail, 
should raise the specter of subsolidus chemical equilibration. 
Another related clue to subsolidus equilibration is that differ- 
ent minerals may imply different magma compositions, un- 
related to order of crystallization (e.g., Fig. 2). 

Thus, the details of mineral chemistry and the appropriate- 
ness of inferred magma compositions may provide important 
clues to whether an igneous rock has experienced significant 
subsolidus equilibration. If there is evidence of significant 
subsolidus equilibration, calculation of magma compositions 
via /) mmcra,/ma g ma involves great peril: the inference of bizarre, 
fractionated magmas from rocks that formed from common, 
simple magmas. 

Note Added in Proof: REE abundances in mineral separates 
from the equilibrated eucrites Stannem and EET79002 con- 
form to the predictions of section 2.1, as do their “parent 
magmas” calculated by D mi " era!/basah (Phinney et al., 1993). 
Phinney W. C., Lindstrom D. J., Mittlefehldt D. W., and Mar- 
tinez R. R. ( 1993) Post-igneous redistribution of components 
in eucrites. Lunar Planet. Sci. XXIV, 1137-1138 (abstr.). 
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